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An  exTjerlnental  study  was  oofPled  out  toJ>lnvestlgatec^— v 
the  actual* Influence  of  the  boundary  layer  on  heat  transfer  ] 

losses  in  typical  expanders  utilized  in  lorr  teanerature  _ I 

systems*--. >fore  specifically,  to  investigate ^e  effecf  of 
the  boundary  layer  on  the  time  denendence  of  the  rate  of 
heat  tr^sfer  through  'tfee^ylinder  walls.  The  experimental 
expanderXtttildbced-  was  a  closed  cylinder  apparatus  with  a 
fixed  mass  charge  of  helium  gas.  Cyclic  pressure  and  volume 
data  points  were  generated  during  operation  of  the  expander. 
The  resulting  experimental  pressure  vs.  volume  data  trace 
was  stored  digitally  and  presented  graphically  using  a 
Nlcolet  storage  oscilloscope.  ^The  CMS  computy  system  was 
used  to  further  reduce  and  interpret  the  data^-T  f  ^ 


The  results^oncluded  that  the  temperature  of  the 
boimdary  layer  gas  had  the  controlling  effect  over  the 
cylinder  bulk  gas  temperature  in  determining  the  actual  rate 
of  heat  transfer  through  the  cylinder  walls.  The  apparent 
zero  rate  of  heat,  transfer  point  occurred  when  the  bulk  gas 
temperature  ecualed  the  cylii.der  w^l_te^erature.-r>  Because 
of  the  strong  effect  that  the  ^ouhdaryT^^  had  on  the 
heat  transfer  in  the  cylinder ,S^e  actual  zero  rats  of  heat 
transfer  point  occurred  when  the  boundary  layer  gas  temper¬ 
ature  equalled  the  cylinder  wall  temperature.  The  boundary 
layer  gasj:emperature  was  determined  responsible  for  up  to  a 
64  degTM^hase  difference  in  the  times  of  the  actual  and 
apparent  zero  rate  of  heat  transfer  points  in  the  compression 
expansion  strokes  in  the  cylinder.^ This  phase  difference 
was  evaluated  to  be  leading  in  both  cases.  Therefore  the 
actual  time  of  zero  rate  of  heat  transfer  preceded  the  ap¬ 
parent  time  of  zero  rate  of  heat  transfer  by  a  significant 
margin  in  each  piston  stroke.  ^  \ 
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ABSTRACT 

An  experimental  study  was  caurrled  out  to  investigate 
the  actual  Influence  of  the  boundary  layer  on  heat  transfer 
losses  in  typical  expanders  utilized  in  low  temperature 
systems*  More  specifically,  to  investigate  the* effect  of 
the  boundary  layer  on  the  time  dependence  of  the  rate  of 
heat  transfer  through  the  cylinder  walls.  The  experimental 
expander  utilized  was  a  closed  cylinder  apparatus  with  a 
fixed  mass  charge  of  helium  gas.  Cyclic  pressure  and  volume 
data  points  were  generated  during  operation  of  the  expander. 
The  resulting  experimental  pressure  vs.  volume  data  trace 
was  stored  digitally  and  presented  graphically  using  a 
Nlcolet  storage  oscilloscope.  The  CMS  computer  system  was 
used  to  further  reduce  and  interpret  the  data. 

The  results  concluded  that  the  temperature  of  the 
boundary  layer  gas  had  the  controlling  effect  over  the 
cylinder  bulk  gas  temperature  in  determining  the  actual  rate 
of  heat  transfer  through  the  cylinder  walls.  The  apparent 
zero  rate  of  heat  transfer  point  occurred  when  the  bulk  gas 
temperature  equaled  the  cylinder  wall  temperature.  Because 
of  the  strong  effect  that  the  boundary  layer  had  on  the 
heat  transfer  in  the  cylinder,  the  actual  zero  rate  of  heat 
transfer  point  occurred  when  the  boundary  layer  gas  temper* 
ature  equaled  the  cylinder  wall  temperature.  The  boundary 
layer  gas  temperature  was  determined  responsible  for  up  to  a 
64  degree  phase  difference  in  the  times  of  the  actual  and 
apparent  zero  rate  of  heat  transfer  points  in  the  compression 
and  expansion  strokes  in  the  cylinder.  This  phase  difference 
was  evaluated  to  be  leading  In  both  cases.  Therefore  the 
actual  time  of  zero  rate  of  heat  transfer  preceded  the  ap¬ 
parent  time  of  zero  rate  of  heat  transfer  by  a  significant 
margin  In  each  piston  stroke. 
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D^TNTTTQN  Q?  SYMBOLS  USZO 


S-ynbol 

Definition 

Units 

A 

heat  transfer  area 

in2 

S 

specific  heat  at  constant 
*  pressure 

3tu/lbjjj-°S 

specific  heat  at  constant 
volume 

3tu/lh  -°R 

31 

h 

surface  coefficient  of  heat 
transfer 

3tu/hr-ft^- 

a 

system  mass 

m 

P 

cylinder  pressure 

Ib^/in^ 

f’o 

reference  cylinder  pressure 

Ib^/in^ 

p 

conp 

compression  stroke  pressure 
at  reference  volume 

lbj./in^ 

espansion  stroke  pressure 
at  reference  volume 

Ib^/in^ 

^mep 

mean  pressure  over  the  com¬ 
pression  stroke  (eqn  21) 

lb  £./in^ 

PSTR 

non-diaensionalized  pressure 
(eqn  22) 

— 

q 

rats  of  heat  transfer 

in-lb sec 

Q 

quantity  of  heat  transferred 

in-lb ^ 

S 

specific  gas  constant 

ft-lb/lb^- 

RPM 

piston  revolutions  per  minute 

rev/nin 

S 

entropy  of  cylinder  hulk  gas 

in-lb 

SSTS  non-dlaenslonalized  entropy 

(eqn  24) 

t  tiae 


sec 


Definition 

Units 

^bulk  gas 

cylinder  bulk  gas  temperature 

reference  temperature 

°R 

T 

boundary 

layer 

boundary  layer  gas  temper¬ 
ature 

^wall 

cylinder  Trail  temperature 

‘’r 

TSTR 

non-dimensionalized  temper¬ 
ature  (eqn  23) 

— 

U 

internal  energy 

in-lb  ^ 

V 

cylinder  volume 

in^ 

reference  cylinder  volume 

in3 

VBDC 

cylinder  volume  at  piston 
bottom  dead  center. position 

in^ 

^/TDC 

cylinder  volume  at  piston 
top  dead  center  position 

ia-5 

VSTR 

non-dimensionalized  cylinder 
volume  (eqn  20) 

-- 

7 

total  cylinder  vrork 

in-lb  ^ 

7 

comp 

7ork  done  during  compression 
stroke 

in-lb 

7 

exp 

work  done  during  expansion 
stroke 

in- lb  .e- 

ratio  of  specific  heats 

(Cj/O^) 

— 

CHAPTER  I 
IHTPQPgCTTQN 


The  reciprocating  piston-cylinder  system  is  a  common 
and  relatively  simple  device  used  in  many  different  applica¬ 
tions  in  today* s  mechanical  world*  Each  day  some  type  of 
reciprocating  piston-cylinder  system  impacts  on  our  lives 
whether  it  be  through  a  reciprocating  automotive  engine  or 
through  some  other  type  of  reciprocating  pumping  device  that 
may  provide  another  personal  service.  An  additional  example 
of  this  type  of  device  is  a  reciprocating  cryogenic  expander, 
^ose  purpose  is  to  produce  a  reduction  in  the  temperature 
of  the  gas  passing  through  it. 

In  all  of  these  types  of  systems  work  is  produced 
through  the  expansion  of  a  fluid  against  a  piston.  For  an 
ideal  expansion  engine  the  expansion  stage  should  be  isen- 
tropic  and  the  work  per  unit  mass  should  be  expressed  as: 

W  =  Cp  (T^  -  T^)  (1) 

However,  because  this  fluid  expansion  Is  not  an  ideal  pro¬ 
cess,  losses  and  inefficiencies  are  evident.  Some  examples 
of  typical  loss  mechanisms  usually  associated  with  through 
flow  expanders  are: 

1*  clearance  losses  -  Because  of  the  clearance 
volume  remaining  when  the  piston  is  at  top  dead  center,  the 
cylinder  cannot  fully  exhaust  the  entire  old  charge.  There- 
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fore  the  fresh  charge  entering  the  piston  is  further  limited 
in  the  actual  volume  that  it  can  ideally  occupy.  The  clear¬ 
ance  losses  then  define  these  variations  from  the  ideal 
case.  The  actuad.  factors  responsible  for  the  clearance  loss 
include  much  more  complex  processes  such  as  blo^r  in-blovr  out 
loss  effects  and  the  intricate  influences  of  the  mass  flo^r 
rate.  These  processes  are  extremely  involved  and  their 
evaluation  here  is  not  required,  so  then  just  calling  atten¬ 
tion  to  their  existence  will  suffice, 

2,  valve  losses  -  Because  the  inlet  and  exhaust 
valves  are  responsible  for  a  resistance  to  the  fluid  flow, 
there  will  be  a  pressure  drop  across  the  valves  that  will 
account  for  some  loss, 

3*  piston  ring  losses  -  There  will  be  some  loss 
associated  with  the  leakage  past  the  piston  rings.  Also 
any  friction  between  the  piston  rings  and  cylinder  walls 
will  generate  some  heat  that  will  be  transferred  to  the 
fluid  and  add  to  the  loss. 

The  above  examples  are  certainly  not  the  only  loss  mech¬ 
anisms  evident  in  a  piston-cylinder  arrangement,  but  they 
are  usually  considered  the  prominent  factors. 

This  experimental  investigation  will  consider  the  loss 
mechanisms  evident  in  a  closed  cylinder.  No  inlet  or  ex¬ 
haust  processes  will  be  evaluated.  The  mass  of  the  working 
fluid  will  remain  constant  throughout  the  compression  expan¬ 
sion  cycle  being  examined.  Therefore  the  causes  of  the 
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losses  evident  in  the  cylinder  will  simplify,  as  the  heat 
transfer  loss  will  he  the  only  prominent  loss  mechanism. 

This  will  then  permit  the  determination  of  the  actual  effect 
that  the  heat  transfer  loss  has  on  the  system. 

Calculation  of  Loss 

To  evaluate  the  thermodynamic  processes  occurring  in 
the  closed  piston-cylinder  the  first  law  of  thermodynamic.- 3 
will  be  utilized.  In  its  most  basic  form  the  first  law 
states  that 


dQ  =  dU  H*  dV/.  (2) 

If  the  pressure  of  the  gas  is  spatially  uniform  at  each 
state  of  the  cycle  and  the  work  done  by  the  gas  is  re¬ 
stricted  to  the  displacement  of  the  piston,  then  the  work 
term  can  be  defined  as 

dW  =  /  PdV,  •  (3) 

For  an  ideal  gas  the  change  in  internal  energy  of  the  gas  is 

dU  s  mC  dT,.  ,,  (k) 

V  -bulk  gas, 

The  first  law  can  then  be  e:cpressed  as 

JpdV  =  dQ  -  mC^  ^^bulk  gas. 

If  the  process  is  assxnsed  to  be  ideal  and  adiabatic, 
then  the  heat  transfer  in  the  system  will  be  zero  (dO  =  0) 
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and  the  first  Ian  nlll  then  simplify  to 

fpdV  3  mC  dTv  (6) 

When  dT  is  eliminated  by  substitution  of  the  ideal  gas 
equation  of  state,  equation  (5)  integrates  to 

y 

PV  3  constant.  (7) 

The  construction  of  the  trace  for  this  process  on  a  pres¬ 
sure  vs.  volume  diagram  will  yield  a  plot  in  which  the  com¬ 
pression  and  expansion  processes  share  the  same  line  and 
therefore  describe  no  enclosed  area  (Figure  1).  There  will 
be  zero  net  work  done  on  the  working  fluid. 

If  the  process  is  not  ideal,  the  compression-expansion 
trace  will  not  follow  the  PV  3  constant  path.  Plotting 
the  non-ideal  expansion  trace  on  a  pressure  vs.  volume 

4 

diagram  will  yield  a  closed  plot  circumscribing  some  en¬ 
closed  area  (Figure  2).  This  enclosed  area  represents  the 
net  work  done  on  the  working  fluid  during  the  cycle  and  also 
represents  the  total  losses  that  the  system  experiences 
during  the  cycle. 

Therefore  the  total  loss  per  cycle  in  the  piston- 
cylinder  arrangement  can  be  expressed  as  the  integral  of 
pressure  times  a  differential  volume  around  the  entire  cycle 

^PdV  3  loss.  (8) 

This  loss  is  specifically  associated  with  just  the  bulk  gas 
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In  the  cylinder*  The  total  loss  per  cycle  can  also  be  ex¬ 
pressed  in  a  different  form  by  considering  the  loss  as  the 
integral  of  the  rate  of  heat  transfer  from  the  bulk  gas 
through  the  wall  around  the  entire  cycle 

dt  s  loss*  (9) 

Boundary  Laver  Effects 

Adjacent  to  the  inside  of  the  cylinder  wall  a  small 
volume  of  gas  with  finite  mass  will  permanently  exist 
through  which  the  heat  flux  must  travel  in  order  to  reach 
the  cylinder  wall*  As  the  piston  moves  in  the  cylinder^  the 
depth  of  this  thin  layer  of  gas  will  vary  constantly  due  to 
well  defined  boundary  layer  effects*  If  the  cyclic  heat 
flow  entering  and  exiting  the  boundary  layer  is  considered, 
the  importance  of  the  boundary  layer  effect  can  be  illus¬ 
trated  (Figure  3)* 

For  an  entire  cycle  the  total  loss  from  the  bulk  gas  to 
the  wall  dt)  should  be  equal  to  the  loss  from  the  bulk 
gas  to  the  boundary  layer  (^q^  dt)  and  should  also  be  equiv¬ 
alent  to  the  loss  from  the  boundary  layer  to  the  wall 
dt )  * 

^q  dt  a  ^q^  dt  =^92 

Because  the  heat  capacitance  of  the  cylinder  wall  is  very 
great  when  compared  to  the  heat  capacitance  of  the  gas, 
the  wall  temperature  can  be  assumed  to  be  constant  when 


SCHEMATIC  OF  BASIC  HEAT  FLOW  WITH  BOUNDAEY  LAYER 

FIGURE  3 
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considering  a  single  cycle. 


If  the  small  volume  of  gas  comprising  the  boundary 
layer  had  no  heat  capacitance  then  at  any  instant  during  the 
cycle  the  rate  of  heat  transfer  from  the  bulk  gas  to  the 
boundary  layer  would  equal  the  rate  of  heat  transfer  from 
the  boundary  layer  to  the  cylinder  wall. 


/^instantaneous,  zero  heat  \ 
q  =  H-|  =  42  I  capacitance  in  boundary  layery  ^  ' 

However,  the  presence  of  this  boundary  layer  would  suggest 
that  It  would  hold  some  finite  value  of  heat  capacitance 
between  the  temperature  extremes  of  the  bulk  gas  and  the 
cylinder  wall.  Assuming  that  the  boundary  layer  does  sup¬ 
port  some  heat  capacitance,  then  as  the  piston  moves  through 
the  cylinder  the  actual  rate  of  heat  transfer  from  the  bulk 
gas  to  the  boundary  layer  should  instantaneously  differ  from 
the  rate  of  heat  transfer  from  the  boundary  layer  to  the 
cylinder  wall. 

/  /  /instantaneous,  finite  heat  "N 

H  ^  q-j  ^  42  I  capacitance  in  boundary  layery  ^ 


The  rate  of  heat  transfer  and  therefore  the  loss  expe¬ 


rienced  from  the  system  depends  heavily  upon  the  temperature 
difference  between  the  working  gas  and  the  cylinder  wall  and 
upon  the  properties  of  the  working  gas. 

q  s  hA  -  T^all^ 
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During  both  the  compression  and  the  expansion  stroke  of  a 
cycle,  the  system  ^11  pass  through  a  point  rrhere  the  gas 
temperature  will  equal  the  cylinder  wall  temperature.  At 
these  two  points  the  rate  of  heat  transfer  will  equal  zero 
as  the  temperature  difference  will  equal  zero  -  *^’7all 

s  0).  Consider  a  system  in  which  the  cylinder  wall  boundary 
layer  has  zero  heat  capacitance.  Here  the  temperature  of 
the  gas  in  the  layer  ad;]acent  to  the  wall  will  be  approxi¬ 
mately  equal  to  the  bulk  gas  temperature.  Therefore  at  the 
points  in  the  cycle  where  the  bulk  gas  temperature  equals 
the  cylinder  wall  temperature,  the  process  will  be  adiabatic 
and  the  rate  of  heat  transfer  will  equal  zero. 

However,  if  the  system  is  assumed  to  have  finite  heat 
capacitance  in  the  thermal  boundary  layer,  the  temperature 
of  the  gas  in  the  layer  adjacent  to  the  wall  could  differ 
significantly  from  the  bulk  gas  temperature.  In  this  system 
the  points  of  zero  rate  of  heat  transfer  will  not  occur 

where  3  "^wall,  ''rtJ.l  occur  at  the  points  where 

T  —  T 

boundary  layer  ”  wall. 

To  illustrate  this  difference,  if  the  compression  ex¬ 
pansion  process  is  plotted  on  a  pressure  vs.  volume  diagram, 
the  points  of  zero  rate  of  heat  transfer  will  occur  where 
the  process  is  adiabatic  at  that  point.  Previously  in  this 
section  it  was  shown  that  the  process  will  be  adiabatic  at 
the  point  where  the  trace  is  tangent  to  a  line  where  PV^  s 
constant.  Figure  4  illustrates  a  pressure  vs.  volume 
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diagram  for  a  system  where  the  points  of  taingency  to  the 
adiabatic  lines  coincide  with  the  points  where  the  bulk  gas 
temperature  equals  the  cylinder  wall  temperature. 

If  a  system  with  finite'  heat  capacitance  In  the 
boundary  layer  Is  considered,  the  points  at  which  the  bulk 
gas  temperature  equal  the  cylinder  wall  temperature  will 
not  coincide  with  the  points  at  which  the  system  trace  Is 
tangent  to  the  adiabatic  lines.  Figure  3  Illustrates  this 
difference.  In  both  the  compression  and  expansion  stroke, 
a  phase  difference  Is  Identified  that  Is  the  difference  (in 
time  or  degrees  of  rotation)  between  the  points  of  actual 
zero  rate  of  heat  transfer  In  the  system  (points  where  the 
process  is  perfectly  adiabatic  where  T^oundary  layer  = 

^wall^  and  the  points  of  apparent  zero  rate  of  heat  transfer 
In  the  system  (points  where 

Therefore  the  effect  of  the  boundary  layer  on  the  heat 
transfer  In  the  cylinder  must  be  determined  In  order  to 
evaluate  the  Impact  of  the  boundary  layer  on  the  total 
system  loss. 

Other  methods  are  certainly  available  that  would  facil¬ 
itate  calculation  of  the  net  cyclic  loss  and  any  phase  dif¬ 
ference  experienced  In  the  points  of  zero  rate  of  heat 
transfer.  One  of  the  most  direct  and  accurate  methods  re¬ 
quires  the  plotting  of  a  temperature  vs.  entropy  diagram 
for  the  cycle.  Figure  6  Is  an  example  of  such  a  trace.  The 
area  enclosed  by  the  curve  represents  the  net  work  done  on 
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PRESSURE  VS,  VOLUME  DIAGRAM  OF  SYSTEM 
',7ITH  ZERO  BOUNDARY  LAYER  HEAT  CAPACITANCE 
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the  working  fluid  which  is  also  equivalent  to  the  total  loss 
experienced  during  the  cycle.  For  any  single  cycle  the 
following  statement  should  be  true: 

^PdV  s^TdS  =  loss.  (14) 

On  this  temperature  vs.  entropy  diagram*  the  points  on  the 
trace  that  are  tangent  to  a  vertical  line  signify  points 
where  the  change  in  entropy  is  equal  to  zero.  At  these 
points  the  process  is  perfectly  adiabatic  auid  signify  points 
of  zero  rats  of  heat  transfer.  The  temperature  of  the  gas 
adjacent  to  the  cylinder  wall  at  these  points  should  be 
equal  to  the  temperature  of  the  cylinder  walls.  If  these 
points  of  actual  zero  rate  of  heat  transfer  differ  from  the 
points  on  the  graph  of  apparent  zero  rate  of  heat  transfer 
(where  *  "^wall^*  difference  exists. 

In  Figure  6  a  phase  difference  is  evident  in  both  the  com¬ 
pression  and  expansion  strokes. 

The  existence  of  this  phase  difference  must  be  proven 
and  then  evaluated  in  order  to  provide  information  necessary 
to  quantify  the  impact  of  the  thermal  boundary  layer  effects 
on  the  total  system  loss. 


CHAPTER  I I 

DESCRIPTTQW  OF  TEST  EOniPMEITT  AND  TESTING  PROCEDURE 

II. A.  Experimental  Expander 

The  apparatus  used  to  investigate  the  theoretical  pro¬ 
cesses  of  the  expander  was  composed  primarily  of  the  motor 
and  cylinder  body  of  a  modified  air  compressor.  The  active 
cylinder  was  bolted  directly  on  top  of  the  original  com¬ 
pressor  cylinder.  This  cylinder  included  the  gas  charging 
port,  temperature  and  pressure  measuring  instruments  and  had 
the  additional  capability  of  water  cooling  the  cylinder,  A 
diagram  of  this  apparatus  is  shown  in  Figure  7,  The  piston 
utilized  to  compress  and  expand  the  gas  was  housed  in  the 
active  cylinder.  The  active  cylinder  bore  was  2,00  inches 
and  the  piston  stroke  was  3*00  inches.  During  operation  the 
piston  remained  within  the  cylinder,  while  its  maximum  up¬ 
ward  travel  at  the  top  dead  center  position  moved  the  piston 
to  the  upper  lip  of  the  cylinder. 

The  cylinder  extension  housed  a  stationary  piston  that 
formed  the  head  of  the  cylinder  volume,  A  threaded  bolt 
extending  through  the  top  of  the  extension  cylinder  and  into 
the  stationary  piston  provided  the  additional  capability  of 
changing  the  volume  ratio.  By  raising  or  lowering  the  sta¬ 
tionary  piston  through  its  four  inch  range  with  this 
threaded  bolt,  the  cylinder  provided  volume  ratios  ranging 
from  1,5  to  2,4« 
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A  strain  gage  pressure  transducer  vras  utilized  to 
measure  cylinder  pressure*  The  pressure  measurement  is 


determined  in  this  type  of  transducer  by  the  fact  that  the 
resistance  of  the  transducer  strain  gages  on  the  pressure 
sensitive  surface  changes  in  proportion  to  the  pressure  on 
the  surface*  The  transducer  output  is  then  proportional  to 
the  cylinder  pressure*  The  transducer  mounted  in  the  cylin¬ 
der  head  vras  of  the  range  0  to  100  psi  and  rras  calibrated 
at  4*342  psi/millivolt  for  a  6.0  volt  '.nput*  In  order  to 
calculate  the  cylinder  pressure,  the  transducer  output  at  a 
knovm  pressure  had  to  be  determined.  Therefore  prior  to 
recording  any  experimental  data,  the  transducer  input  at 
atmospheric  pressure  rras  recorded  to  provide  the  required 
pressure  reference*  To  make  this  step  easier  a  balar.o'r.  box 
7/as  installed  rrith  the  pressure  transducer  ':/hich  permitted 
the  setting  of  a  known  pressure  directly  on  the  oscilloscope 

II. C.  Volume  Measurement 

The  method  utilized  to  measure  the  cylinder  volume 
consisted  of  a  circular  cam  mounted  eccentrically  on  the 
compressor  crankshaft*  Because  each  cycle  of  the  compressor 
piston  corresponded  to  one  revolution  of  the  compressor 
crankshaft,  the  angle  of  a  given  spot  on  the  circumference 
of  the  shaft  was  directly  proportional  to  the  position  of 
the  compressor  piston  and  therefore  was  also  proportional  to 


the  cylinder  volume*  A  follovrer  was  mounted  on  the  outer 
diameter  of  the  circular  cam  so  that  the  displacement  of  the 
follower  was  proportional  to  the  displacement  of  the  piston* 
This  system  is  displayed  in  Figure  8* 

An  HP  7DCT  displacement  transducer  was  connected  to  the 
follower  and  provided  a  D*C.  output  that  was  proportional  to 
the  displacement  of  the  follower.  The  calibration  of  the 
transducer  provided  the  maximum  output  signal  when  the  cylin¬ 
der  volume  was  maximum  at  the  piston  bottom  dead  center 
position*  3y  knowing  the  maximum  and  minimum  cylinder  volume 
for  the  experimental  run,  the  intermediate  values  for  cylin¬ 
der  volume  could  easily  be  calculated  by  using  the  transducer 
output* 

Cylinder  ’Tall  Tenneraturs  Measurement 
Four  thermocouples  were  utilized  to  measure  the  average 
cylinder  wall  temperature  during  testing*  Three  thermocou¬ 
ples  monitored  the  temperature  at  different  locations  on  the 
outside  of  the  cylinder  wall*  The  remaining  thermocouple 
was  placed  on  the  outside  of  the  cylinder  cooling  water 
Jacket,  The  output  of  the  thermocouples  was  displayed  using 
an  Auto data  Nine  Digital  data  logger  that  provided  a  contin¬ 
uous  readout  in  degrees  Centigrade* 


Digital  Storage  Oscilloscope 
The  output  of  the  pressure  and  volume  transducers  was 
provided  as  input  to  a  Nicolet  Digital  Oscilloscope  Model 
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206 •  The  signal  inforaation  provided  to  the  scope  \7as  con¬ 
verted  to  digital  form,  stored  in  a  buffer  menory  and  then 
vras  displayed  graphically.  Digital  information  for  each 
point  remained  available  in  storage.  The  Nicolet  scope  has 
the  capability  of  storing  up  to  4096  data  points.  Because 
both  pressure  and  volume  were  being  monitored,  each  storage 
cycle  consisted  of  2048  pressure  points  and  2048  volume 
points.  The  scope  had  the  capability  of  graphic  display  of 
either  the  pressure  or  volume  trace  vs.  time  or  an  X-Y  axis 
plot  of  pressure  vs.  volume.  The  tine  interval  between  data 
points  could  be  chosen  on  the  scope  to  ensure  that  an  entire 
cycle  of  data  would  be  available  in  a  single  storage  cell. 


II. F.  Testing  Procedure 

The  following  steps  were  performed  in  this  order  during 
the  data  acquisition  run: 

1.  The  position  of  the  stationary  piston  was  set 
properly. 

2.  The  piston  was  moved  to  bottom  dead  center 
position  and  then  checked  by  assuring  that 
the  output  of  the  volume  transducer  was  at  a 
ma^mum. 

3.  The  cylinder  was  opened  to  the  atmosphere  and 
the  pressure  transducer  output  was  set  on  the 
oscilloscope  and  recorded. 

4.  The  cylinder  was  charged  with  gas  and  the 
pressure  transducer  output  recorded.  At  this 
point  must  ensure  that  the  cylinder  is  not 
charged  to  a  pressure  that  will  over-pressurize 
the  pressure  transducer  during  compression. 

5.  The  compressor  motor  was  turned  on  and  the 
oscilloscope  set  to  automatically  monitor  a 
single  cycle  of  the  cylinder. 
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The  cylinder  wall  temperatures  were  recorded 
Immediately  upon  completion  of  the  sample 
cycle. 

The  compressor  motor  was  secured. 

The  pressure  7s.  volume  data  was  monitored  on 
the  oscilloscope  to  ensure  that  the  trace  was 
satisfactory. 


.  t/r:  mr 


*•  5 
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In  order  to  efficiently  interpret  and  display  the  raw 
data  extracted  from  the  Nicolet  oscilloscope,  the  Institute's 
CMS  computer  system  was  used*  The  data  was  entered  into  the 
computer  files  and  then  was  operated  on  to  provide  various 
thermodynamic  parameters  at  each  point*  Using  this  infor¬ 
mation  thermodynamic  traces  for  the  cycle  were  plotted* 
Finally  the  losses  for  the  cycle  were  calculated  by  numer¬ 
ical  integration  of  the  curves* 

TIT. A.  Computer  Inruts 

The  digital  pressure  vs*  volume  data  provided  by  the 
Nicolet  oscilloscope  was  input  into  a  data  file  that  was 
accessible  to  the  main  program*  Other  initial  inputs  re¬ 
quired  included: 

1*  The  type  of  working  gas  used*- 

£•  The  specific  gas  constant  for  the  gas  (H). 

3*  The  specific  heat  at  constant  pressure  for 

th.  (Cp). 

The  specific  heat  at  constant  volume  for  the 
gas  (C^)* 

5*  The  cylinder  volume  at  the  top  dead  center 
position  (VTDO* 

6*  The  cylinder  volume  at  the  bottom  dead  center 
position  (VBDC). 

7#  The  total  number  of  data  points  (?,  V  pairs) 
entered  into  the  computer  data  file. 

d*  The  time  interval  between  data  points. 
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In  order  to  have  the  capability  of  accurately  coraparing 
data  at  varying  volune  ratios,  compressor  speeds,  and  for 
different  gases,  computations  performed  on  the  data  were  done 
using  common  reference  points*  Therefore  a  standard  method 
'nas  required  to  determine  each  of  the  required  temperature, 
volume  and  pressure  reference  points* 

Cylinder  bulk  gas  temperature  was  the  first  reference 
point  determined  and  was  set  at  100°?.  This  temperature  'vas 
chosen  as  it  was  a  typical  averaige  cyclic  temperature  experi¬ 
enced  in  previous  experimentation  done  with  the  same  e:raeri- 
mental  cylinder  apparatus*  ‘^hen  expressed  on  the  absolute 
scale  the  reference  temperature  was  then  set  at  the  following 
value: 

Tg  «  559.67°  R 

The  method  used  to  determine  the  reference  volume  con¬ 
sisted  of  calculating  the  geometric  mean  using  the  volume  of 
the  cylinder  at  the  top  dead  center  and  bottom  dead  center 
positions*  The  actual  formula  for  the  determination  of 
reference  volume  was 
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Because  the  cylinder  volume  v/as  expressed  in  units  of  cubic 
inches,  the  reference  volume  vas  also  in  cubic  inches. 

The  determination  of  the  reference  pressure  point 
required  the  prior  determination  of  the  reference  volume 
point.  Using  the  reference  volume  point,  the  pressures 
corresponding  to  this  volume  for  both  the  compression  and 
expansion  strokes  were  found.  The  reference  pressure  was 
then  determined  to  be  the  average  of  the  compression  and 
expansion  pressure  values  at  the  reference  volume. 


p  +  p 
com-Q  e:co 


(16) 


Reference  pressure  was  a:qjressed  in  units  of  pounds  force 
per  square  inch. 


In  order  to  calculate  the  bulk  gas  temperature  and 
entropy  at  each  data  point  the  system  mass  had  to  be  deter*- 
mined.  This  calculation  was  performed  using  the  ideal  gas 
equation  of  state  and  the  reference  temperature,  volume,  and 
pressure  values. 


m  -  0  0 


(17) 


The  system  mass  was  then  expressed  in  units  of  pounds  mass, 


Temwerature  and  Sntro-nv 


The  bulk  gas  temperature  at  each  data  point  was  then 


determined  utilizing  the  ideal  gas  equation  of  state  and 


other  previously  determined  parameters. 


bJL 

R  m 
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(18) 


The  temperature  values  were  expressed  on  the  absolute  scale 
in  degrees  Sankine. 

The  entropy  value  for  each  data  point  was  calculated 
referenced  to  the  previously  determined  values  of  reference 
volume  and  pressure  by  the  relation 

-  Sq  s  AS  s  mC^  In 

The  entropy  values  were  expressed  in  units  of  inch-pounds 
force  per  degree  Rankine. 

III.S.  Calculation  of  System  Losses 

In  order  to  calculate  the  system  losses  during  the 
cycle,  cyclic  traces  of  pressure  vs.  volume  and  temperature 
vs«  entropy  were  plotted.  The  area  enclosed  within  either 
of  these  closed  curves  should  equal  the  total  system  loss. 

^PdV  a^TdS  =  loss  (lif) 

The  integration  was  performed  numerically  using  the  trap- 

/ 

ezoid  rule  at  each  successive  data  point.  Total  area  under 
the  expansion  curve  was  considered  negative.  The  total 
system  loss  was  then  described  in  the  units  of  negative 


+  mC^  In 


^0 


(19) 
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inch-pounds  force. 

An  attempt  vras  made  to  calculate  the  area  using  a  much 
more  involved  and  typically  more  accurate  five  point  Gaus¬ 
sian  quadrature  integration  method.  Hor;ever,  because  of 
the  scatter  of  the  digital  data  points  near  the  ends  of  the 
curves,  this  method  proved  no  more  accurate  than  the  much 
simpler  trapezoid  integration  method. 

III.F.  Non-Dimensional  Parameter  Calculations 

In  order  to  provide  the  additional  capability  of  com¬ 
paring  the  computer  output  to  the  output  data  of  any  other 
similar  apparatus,  notxTithstanding  differences  in  the  units 
involved,  non-dimensional  parameters  vrere  introduced.  The 
computer  routine  then  also  included  the  calculations  re¬ 
quired  to  produce  the  output  of  non-dimensional  values  of 
pressure,  volume,  temperature  and  entropy.  Plots  were  also 
made  of  non-dimensional  pressure  vs.  volume  and  non-dimen¬ 
sional  temperature  vs.  entropy  in  order  to  determine  the 
non-dimensional  loss  values. 

III.F.1.  Non-Dimensional  Volume 

The  volume  data  points  were  non-dimensionalized  by 
dividing  each  data  point  by  the  total  displaced  volume  of 
the  cylinder. 


VSTB^  s 


C'/Bic-'/TDC} 


^20) 
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The  cylinder  pressure  data  points  were  non-dimension- 
alized  by  dividing  each  data  point  by  a  mean  pressure  cal¬ 
culated  from  the  compression  stroke.  This  mean  pressure  was 
defined  as  the  total  work  added  during  the  compression 
stroke  divided  by  the  total  displaced  volume  of  the  cylinder 

7/ 

o  _  como  *  ^  o  i  > 

^mep  “  (^/Sdc-v^:dc) 

Each  pressure  data  point  was  then  divided  by  this  mean  pres¬ 
sure  during  compression  to  find  the  non-dimensional  pressure 
values. 


PSTH^ 


(22) 


Non-Dimensional  Temperature 

The  bulk  gas  temperature  data  points  v;ere  non-dimension 
alized  by  dividing  each  data  point  by  the  reference  temper¬ 
ature. 


TSTH.  rr  si  (23) 

0 

III.F.h.  Non-Dimensional__EntroT3y 

The  entropy  values  were  non-dimensionalized  by  dividing 
the  calculated  entropy  values  by  the  mean  pressure  during 
compression  and  the  displaced  volume  of  the  cylinder  and 
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then  multiplying  by  the  reference  temperature. 


Non-dimensional  values  of  total  system  loss  Trere  cal¬ 
culated  in  the  same  manner  as  the  dimensional  loss  values 
were  determined.  Plots  were  made  of  non-dimensional  pres¬ 
sure  vs.  volume  and  non-dimensional  temperature  vs.  entropy 
and  the  enclosed  areas  were  calculated  by  numerical  inte¬ 
gration  utilizing  the  trapezoid  method. 


III.G.  Corn-outer  Outnut 

A  listing  of  the  Fortran  computer  program  is  displayed 
in  Appendix  A.  A  listing  of  the  test  data  file  used  in  the 
computations  is  presented  in  Appendix  B. 

The  output  of  the  computer  routine  includes  the  fol¬ 
lowing  data: 

1.  Pressure,  volume,  bulk  gas  temperature  and 
entropy  at  each  data  point. 

2.  Non-dimensional  pressure,  volume,  bulk  gas 
temperature  and  entropy  at  each  data  point. 

3«  Piston  revolutions  per  minute. 

4*  Time  interval  between  data  points. 

5*  Reference  temperature  volume  and  pressure 
values. 

6.  Average  cylinder  wall  temperature, 

7*  System  mass. 


8.  Pressure  vs,  volume  plot. 

9.  Temperature  vs.  volume  plot. 

10.  Temperature  vs,  entropy  plot. 

11.  Non-dimensional  pressure  vs.  non-dimensional 
volume  plot, 

12.  Non-dimensional  temperature  vs.  non-dimen¬ 
sional  volume  plot,* 

13.  Non-dimensional  temperature  vs,  non-dimen¬ 
sional  entropy  plot, 

14.  Calculations  of  dimensional  and  non-dimen¬ 
sional  system  loss. 


An  example  of  the  computer  output  is  displayed  in  Appendices 
C  and  D, 


CHAPTER  IV 
RESULTS 


I  The  results  of  the  single  data  set  used  are  presented 

*  in  the  computer  outnut  in  Apnendices  C  and  D.  Helium  iras 

I  .  ... 

f  .  utilized  in  the  cylinder  as  the  TOrking  gas  for  the  exuer- 

■ 

iment* 

The  final  pressure  vs,  volume  and  non-dimensional 
pressure  vs,  non-dimensional  volume  plots  resulted  in 
extremely  smooth  traces,  Each  depicted  what  could  easily 
be  identified  as  a  conventional  pressure  vs,  volume  trace 
for  this  type  of  apparatus.  Basically  the  same  comment 
could  be  made  about  both  the  temperature  vs,  volume  and 

j 

I  non-dimensional  temperature  vs,  non-dimensional  volume 

plots.  However,  a  greater  degree  of  scatter  seemed  to  be 
evident  in  the  temperature  vs,  volume  plots  than  in  the  . 
pressure  vs,  volume  plots. 

In  the  temperature  vs,  entropy  and  non-dimensional 
temperature  vs,  non-dimensional  entropy  plots  the  scatter 
is  much  more  evident.  The  shape  of  the  trace  also  is 
somewhat  different  from  a  shape  that  one  might  expect  the 
plot  to  have  after  first  seeing  the  pressure  vs,  volume  and 
temperature  vs.  volume  traces.  Figure  9  is  a  single  line 
curve  approximation  of  the  TSTR  vs,  SSTR  computer  trace 
from  Appendix  D,  The  figure  shows  the  adiabatic  and  iso¬ 
thermal  limits  of  a  plot  of  this  type.  The  adiabatic  limit 
being  a  curve  where  the  compression  and  expansion  strokes 
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TSTH  VS.  SST!?  SMOOTH  CUHVZ  THACE  jITK  ADIABATIC  AlTD  ISOTHEHI^LAL 
LIMITS  AND  COMPHESSION  AiTD  EXPANSION  STSOKE  PHASE  DIEFEP2D:CES 


FIGTOE  9 


approach  being  perfectly  adiabatic  while  the  isothernal 
limit  approaches  a  curve  where  the  compression  and  erpauision 
strokes  are  perfectly  isothermal.  The  TSTS  vs.  SSTH  plot 
has  been  drawn  on  a  graph  with  auces  of  equal  length  so  that 
the  curve  has  not  been  mistakenly  lengthened  in  either 
direction  to  give  the  false  appeeurance  of  being  more  adia¬ 
batic  or  isothermal.  Nevertheless,  the  trace  seems  to 
approach  the  adiabatic  limit  in  the  high  temperature  portion 
of  the  plot,  while  approaching  the  isothermal  limit  in  the 
low  temperature  portion  of  the  plot.  On  the  compression 
stroke  the  trace  becomes  perfectly  adiabatic  at  a  point 
early  in  the  stroke  at  TST2  s  .365  (T  =  On  the 

expansion  stroke  the  trace  also  becomes  perfectly  adiabatic 
at  a  point  early  in  the  stroke  at  TST2  »  1.15  (T  s  644° H) . 

At  each  of  these  two  points  that  actual  rate  of  heat 
transfer  in  the  system  through  the  cylinder  walls  equals 
zero. 

The  points  of  apparent  zero  rate  of  heat  transfer  occur 
at  T  =  "^wall  ”  539 The  difference  between  these  points 
and  the  points  of  actual  zero  rate  of  heat  transfer  is  the 
phase  difference.  In  both  the  compression  and  expansion 
strokes  the  phase  difference  is  leading  as  the  point  of 
actual  zero  rate  of  heat  transfer  occurs  prior  to  reaching 
the  point  of  apparent  zero  rate  of  heat  transfer.  The 
actual  phase  differences  can  be  computed  to  be: 
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compression  stroke  k3  leading 
expansion  stroke  64®  leading 

The  actual  system  loss  figures  computed  were: 

PdV  loss  a  -24.951  in-lb ^ 

TdS  loss  a  -25.107  in-lb^ 

Non-dimensional  PdV  loss  a  -0.0584 
Non-dimensional  TdS  loss  a  -0.0588 

In  both  the  dimensional  and  non-dimensional  cases,  the 
and  TdS  losses  showed  less  than  one  percent  difference 


CHAPTS5  V 

discussion 


V.A.  Data  Point  Scatter 

The  scatter  evident  in  the  temperature  vs.  volume  plots 
and  particularly  in  the  temperature  vs.  entropy  plots 
deserves  explanation.  In  these  plots  this  scatter  occurs 
primarily  in  the  center  section  of  the  trace,  away  from  the 
aureas  of  minimum  and  maximum  cylinder  volume.  At  the  upper 
and  lower  end  portions  of  all  the  plots,  the  data  points  are 
generally  packed  tightly  together  in  an  orderly  fashion. 

This  is  due  to  the  fact  that  at  the  ends  of  the  compression 
and  expansion  stroke,  the  piston  is  decelerating  rapidly, 
halting  abruptly  and  then  again  rapidly  accelerating  away 
from  its  extreme  position.  This  stopping  and  starting 
motion  requires  a  great  deal  of  time  to  cover  a  small  dis¬ 
tance  in  comparison  to  the  piston' s  much  more  rapid  pace 
at  the  mid-stroke  position  of  its  travel.  Because  the  data 
points  are  taken  at  equal  time  intervals,  many  more  data 
points  are  taken  at  the  extreme  positions  of  the  plots 
providing  the  tighter  packing  and  the  generally  smoother 
appearance  in  those  areas. 

Near  the  mid-stroke  positions  of  the  cycle  a  different 
situation  exists.  Here  the  piston  is  moving  much  more 
rapidly  relative  to  the  velocity  experienced  at  the  extreme 
cylinder  positions.  Again  as  the  data  points  are  taken  at 
equal  time  intervals,  a  greater  distance  should  exist 
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bet'sreen  points  on  each  trace  near  the  mid-stroke  position. 
This  situation  alone  does  not  account  for  the  scatter,  but 
along  ^th  t'^ro  other  possible  causes,  the  existence  of  the 
scatter  can  be  better  explained. 

First,  the  Nicolet  oscilloscope  samples  2043  volume 
and  2948  pressure  data  points  per  scope  storage  cycle  at 
equal  time  intervals  between  each  point.  The  procedure  of 
the  scope  is  to  sample  one  volume  data  point,  wait  the 
appropriate  time  interval  and  then  sample  one  pressure  data 
point,  continuing  in  this  manner  until  the  storage  cycle  is 
full.  The  scope  does  not  have  the  capability  to  sample 
pressure  and  volume  simultaneously.  Therefore  the  resulting 
pressure- volume  data  point  pair  formed  by  matching  a  volume 
data  point  with  the  following  pressure  data  point  does  not 
exactly  coincide  with  that  pressure-volume  data  point  pair 
that  would  occur  if  the  capability  existed  to  sample  each 
pair  simultaneously.  This  fact  brings  a  small  amount  of 
error  into  the  position  of  each  data  point. 

Second,  the  digital  nature  of  the  data  generates  addi¬ 
tional  error.  The  accuracy  of  the  digital  data  output  of 
the  Nicolet  was  limited  to  .05  millivolt.  7ith  the  large 
number  of  data  points  taken,  the  changes  in  output  voltage 
between  many  points  was  less  than  .05  millivolt.  Therefore 
voltage  changes  of  less  than  .05  millivolt  were  rounded  off 
to  the  next  .05  millivolt  increment  creating  some  error  in 
the  data.  This  might  have  not  been  critical  except  for  the 
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fact  that  the  volume  data  was  being  read  on  the  ICO  milli¬ 
volt  scale  on  the  oscilloscope,  and  the  range  of  the  data 
covered  only  3*00  -  17.90  millivolts  (roughly  15/3  of  the 
range  of  the  scale).  The  100  millivolt  scale  vras  the  best 
scale  to  use  as  it  was  the  smallest  scale  available  on  the 
scope  that  would  accept  the  range  of  the  volume  data.  How¬ 
ever,  a  or  5X  signal  amplifier  should  have  been  connected 

in  the  volume  signal  wiring  to  utilize  the  full  range  of  the 
100  millivolt  scale.  This  would  have  alleviated  the  problem 
existing  due  to  the  .05  millivolt  accuracy  of  the  riicolet 
oscilloscope. 

*  These  limitations  of  the  digital  oscilloscope  in  con¬ 
junction  with  the  magnitude  of  the  mid-stroke  piston  velocity 
producing  a  greater  distance  between  data  points  in  this 
area  of  the  plots,  are  surely  a  sufficient  basis  with  which 
to  attempt  to  justify  the  scatter. 

V.B.  Shane  of  Temwerature  vs.  Entrbov  Plots 

The  actual  shape  of  the  temperature  vs.  entropy  plots 
exhibits  a  few  variations  from  what  might  be  expected  after 
seeing  the  other  plots.  In  the  high  temperature  end  of  the 
plot,  the  curve  seems  to  approach  its  adiabatic  limit, 
especially  on  the  initial  portion  of  the  expansion  stroke 
(Figure  10).  The  explanation  for  this  is  fairly  simple. 

At  the  start  of  the  expansion  stroke,  the  temperature  of  the 
bulk  gas  and  the  temperature  of  the  gas  in  the  boundary 
layer  are  both  greater  than  the  cylinder  wall  temperature. 
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However,  the  temperature  of  the  bulk  gas  is  still  signifi¬ 
cantly  greater  than  the  boundary  layer  gas  because  the  lower 
temperature  cylinder  wall  has  a  greater  effect  on  the  temper 
ature  of  the  gas  in  the  thin  boundary  layer.  As  the  expan¬ 
sion  stroke  begins,  the  pressure  drops  and  the  bulk  gas 
temperature  drops  rapidly.  The  temperature  of  the  boundary 
layer  gas  drops  also,  but  it  drops  at  a  much  slower  rate  as 
it  started  the  expansion  stroke  at  a  significantly  lower 
temperature.  Because  the  temperature  difference  between  the 
boundary  layer  gas  and  the  cylinder  7/all  is  not  great  for 
the  first  thirty  degrees  of  the  expansion  stroke,  very  littl 
heat  is  transferred  out  of  the  system  and  the  process  ap¬ 
proaches  its  adiabatic  limit. 

At  the  low  temperature  end  of  the  plot,  near  the  end 
of  the  expansion  stroke  the  curve  seems  to  approach  its 
isothermal  limit  (Figure  10).  The  explanation  for  this  is 
very  similar  to  that  for  the  near  adiabatic  case.  Near  the 
end  of  the  expansion  stroke  the  boundary  layer  gas  temper¬ 
ature  has  fallen  somewhat  lower  than  the  cylinder  wall  tem¬ 
perature.  In  the  last  thirty  degrees  of  the  acpansion  stroke 
the  volume  in  the  cylinder  increases,  but  it  increases  at  a 
relative  rate  much  less  than  the  rate  of  increase  near  the 
start  of  the  expan;^on  stroke.  Therefore  the  change  in  the 
bulk  gas  temperature  will  also  change  at  a  rate  much  less 
than  the  rate  of  change  at  the  start  of  the  expansion  stroke 
Sven  though  the  process  approaches  its  isothermal  limit. 


there  is  still  a  sufficient  temperature  difference  betv/een 
the  boundary  layer  gas  and  C3^1inder  ?rall  to  permit  signifi¬ 
cant  heat  transfer  from  the  system. 

The  temperature  vs.  entropy  plot  also  has  v;hat  seems  to 
be  an  inflection  near  the  mid-stroke  position  on  the  compres¬ 
sion  stroke  (Figure  10).  At  this  inflection  there  is  an 
abrupt  change  in  the  slope  of  the  compression  trace.  Because 
of  its  proximity  to  the  mid-stroke  position  on  the  curve, 
this  inflection  may  be  due  to  the  change  of  piston  accelera¬ 
tion  to  deceleration  at  mid-stroke.  No  such  glaring  inflec¬ 
tion  is  evident  in  the  expansion  curve, 

7.C.  Identification  of  the  Adiabatic  Points  in  the  Cvcle 

The  identification  of  the  tvro  points  in  the  cycle  that 
are  adiabatic  is  extremely  important  in  determining  the 
actual  phase  difference  bettreen  the  heat  transfer  and  the 
bulk  gas  to  '.Tall  temperature  difference.  It  is  e::pected  that 
both  the  compression  and  expansion  strokes  should  each  have 

one  adiabatic  point  vbere  Tboundary  layer  =  ^oall. 
this  should  occur  at  the  points  of  ma:d.mum  and  minimum  values 
of  entropy.  However,  because  of  the  scatter  in  the  tempera¬ 
ture  vs.  entropy  plots  it  was  difficult  to  envision  the 
actual  location  where  the  curve  became  tangent  to  the  ver¬ 
tical.  Therefore  to  identify  the  location  a  rough  curve  was 
passed  through  the  data  points  (Figure  9)  and  the  adiabatic 
points  were  then  evaluated  at  the  locations  ’.There  the  curve 
passed  through  the  vertical. 


To  check  these  adiabatic  locations,  lines  of  ?V  = 
constant  (adiabatics)  ^ers  constructed  on  the  pressure  vs. 
volume  plot  (ITisure  n).  The  locations  on  the  trace  that 
were  tangent  to  an  adiabatic  line  vrould  signify  an  adiabatic 
point.  The  results  of  this  procedure  verified  the  selection 
of  the  adiabatic  points  in  Figure  9. 

V.D.  Effects  of  the  Boundary  Laver 

The  effects  of  the  boundary  layer  on  the  rate  of  heat 
transfer,  or  more  specifically  on  the  losses  in  the  system 
are  most  important.  It  is  obvious  that  it  is  the  temperature 
of  the  boundary  layer  gas  and  not  the  bulk  gas  temperature 
^Thich  is  the  driving  force  on  the  transfer  of  heat  through 
the  cylinder  Trail.  If  the  boundary  layer  gas  temperature  is 
greater  than  the  temperature  of  the  cylinder  ’vall,  the  system 
^11  lose  heat.  If  the  boundary  layer  gas  temperature  is 
less  than  the  temperature  of  the  cylinder  vrall,  the  system 
will  gain  heat. 

Along  the  compression  stroke  and  part  of  the  expansion 
stroke  from  the  minimum  to  maximum  entropy  value  the  boundary 
layer  acts  to  cool  the  system  (Figure  10).  On  this  area  of 
the  curve  the  boundary  layer  gas  temperature  is  basically 
greater  than  the  cylinder  wall  temperature  and  heat  vhLll  be 
transferred  out  of  the  system  through  the  cylinder  walls. 

The  heat  transfer  in  this  portion  of  the  curve  is  driven  by 
the  change  in  volume  and  not  the  change  in  the  bulk  gas  tem¬ 
perature  of  the  system. 
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Along  the  8:cpansion  stroke  froa  the  niniaum  entropy 
value  to  the  masiaum  entropy  value  the  boundary  layer  acts 
to  heat  the  system  (Figure  10),  The  temperature  of  the  bulk 
gas  and  the  boundary  layer  is  falling;  hovvever,  the  boundary 
layer  temperature  is  generally  less  than  the  cylinder  v/all 
temperature  and  therefore  heat  nill  be  transferred  into  the 
system  through  the  cylinder  vralls, 

y.S.  Zxnlanation  of  the  Phase  Difference 

The  phase  difference  that  exists  in  the  zero  rate  of 
heat  transfer  location  is  a  very  important  effect  of  the 
boundary  layer  (Figure  9)*  The  e:q>lanation  of  the  cause  for 
the  phase  difference  should  be  fully  discussed. 

The  bulk  gas  temperature  is  an  important  parameter  in 
this  type  of  system,  Because  the  bulk  gas  temperature  is 
much  easier  to  measure  and  to  understand  than  the  more  dif¬ 
ficult  to  measure  boundary  layer  temperature,  the  bulk  gas 
temperature  is  utilized  much  more  in  this  type  of  system. 
Therefore,  if  the  bulk  gas  temperature  is  utilized,  then 
vrhen  deriving  the  apparent  expression  that  vrould  define  the 
rate  of  heat  transfer  from  the  system,  the  temperature  ex¬ 
tremes  utilized  ’vould  be  the  bulk  gas  temperature  and  the 
cylinder  ’^all  temperature  and  the  expression  would  be: 

q  s  hA  gag  “  '^vfall^* 

For  this  eirjression  the  points  of  zero  rate  of  heat  transfer 
are  the  points  on  the  cycle  where  gag  =  ^-all. 


The  existence  of  a  boundary  layer  ivith  finite  heat 
capacitance  altars  the  former  interpretation  of  the  problem. 
The  temperature  of  the  gas  in  the  layer  next  to  the  cylinder 
trail  trill  not  equal  the  bulk  gas  temperature,  but  trill  equal 
the  boundary  layer  gas  temperature.  Therefore,  the  expres¬ 
sion  defining  the  rate  of  heat  transfer  to  the  cylinder  trails 
trill  become: 

q  3  hA  ^^laoundary  layer  “  'rrall^* 

This  expression  specifies  that  the  points  of  zero  rate  of 
heat  transfer  trould  be  the  points  on  the  cycle  trhere 

‘boundary  layer  =  '^•Jall.  <^>'®  S®®  temperature 

is  not  equal  to  the  boundary  layer  gas  temperature  in  this 

system,  some  difference  trill  exist  in  the  time  of  apparent 
zero  rate  of  heat  transfer  and  actual  zero  rate  of  heat 
transfer.  This  difference  is  the  zero  rate  of  heat  transfer 
phase  difference.  The  phase  difference  in  both  the  com¬ 
pression  and  expansion  strokes  of  this  experimentation  vras 
leading,  or  the  actual  zero  rate  of  heat  transfer  occurred 
in  the  cycle  before  the  apparent  zero  rate  of  heat  transfer 
(Figure  9). 

The  explamation  for  the  phase  difference  e:d.sts  in  the 
examination  of  the  temperature  history  of  the  boundary  layer. 
As  the  compression  stroke  is  started  the  boundary  layer  gas 
temperature  is  slightly  greater  than  the  bulk  gas  temperature. 
This  is  due  to  the  fact  that  the  cylinder  ’rrall  temperature  is 
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greater  than  the  bulk  gas  temperature  and  the  cylinder  Trail 
temperature  has  a  much  greater  effect  on  the  temperature  of 
the  boundaury  layer  gas  rather  than  the  bullc  gas.  As  the  gas 
is  compressed  the  bulk  gas  and  boundary  layer  gas  temperatur 
rise,  but  because  the  boundary  layer  gas  temperature  started 
the  stroke  at  a  higher  temperature  than  the  bulk  gas,  the 
boundary  layer  gas  temperature  passes  through  the  cylinder 
wall  temperature  prior  to  the  bullc  gas  providing  a  leading 
phase  difference. 

As  the  erpansion  stroke  is  started  both  the  bulk  gas 
and  boundary  layer  gas  are  cooled  by  the  drop  in  pressure, 
3ut  because  the  boundary  layer  gas  temperature  was  signifi¬ 
cantly  lox7sr  than  the  bulk  gas  temperature  at  the  beginning 
of  the  expansion  stroke,  the  boundary  layer  gas  temperature 
passes  through  the  cylinder  wall  temperature  prior  to  the 
bulk  gas  providing  a  leading  phase  difference. 


52 


CHAPTER  VI 


CQHCLJSIQHg 

One  of  the  original  objectives  of  the  vrork  included 
the  capability  of  direct  transzaission  of  the  data  from  the 
experimental  expander  to  the  Institute's  CT'!S  computer  system. 
This  capability  required  the  use  of  a  converter  that  could 
properly  interface  vrith  the  Hicolet  storage  oscilloscope  and 
the  CMS  computer  system.  The  Electric  Pov;er  Systems  Labora¬ 
tory  was  assigned  the  task  of  designing  and  assembling  this 
interface  unit.  Due  to  manpower  shortages,  the  Electric 
Power  Systems  Laboratory  was  not  able  to  complete  the  as¬ 
sembly  of  the  interface  unit  prior  to  the  submission  of  this 
paper.  The  unavailability  of  the  interface  unit  required  the 
data  to  be  read  point  by  point  digitally  from  the  Kicolet 
storage  oscilloscope  and  then  manually  entered  into  a  CI'iS 
computer  system  data  file. 

'7ith  the  use  of  an  operable  interface  unit,  dozens  of 
cycle  runs  could  have  bean  performed.  However,  the  lack  of 

the  interface  significantly  limited  the  number  of  cycle  runs 
that  could  be  evaluated.  Further  setbacks  occurred  '"hen 
serious  problems  were  encountered  in  the  digital  and  graphic 
data  presentation  from  the  Hicolet  storage  oscilloscope.  The 
final  result  was  the  availability  of  only  one  good  set  of 
cycle  data  for  analysis. 

Although  the  technique  developed  in  this  paper  to 
analyze  the  data  was  used  on  only  one  set  of  experimental 
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data,  the  results  obtained  proved  that  the  analysis  nethod 
rras  highly  credible*  '.Vhether  one  or  one  hundred  sets  of 
cycle  data  had  been  evaluated,  the  sane  method  ?;ould  have 
been  utilized* 

The  additional  runs  that  could  have  been  performed  ’.vith 
an  operable  interface  unit  would  have  greatly  increased  the 
scope  of  the  experiment*  Parameters  such  as  piston  HPM, 
Heynold's  number,  cylinder  volume  ratio,  cylinder  pressure 
and  type  of  working  gas  could  have  bean  varied.  The  results 
that  would  have  been  presented  with  the  variations  in  the 
above  parameters  would  have  provided  insight  into  the  actual 
heat  transfer  losses  cind  phase  difference  for  a  '.Tide  range 
of  conditions*  This  general  information  could  then  be 
extranolated  for  use  in  the  design  of  other  devices  of  this 
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